function [bmean]=olsmoment(Y,X1)

bmean=(X1'*X1)\(X1'*Y);


